Some spectral properties of the quasar ultraviolet bump 
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Abstract. In the present work the part of the quasar UV-optical bump within 
the wavelength range 1210- 1450 A was studied with the help of composite spectra 
compiled from the samples of SDSS DR7 spectra with the similar spectral index 
a\ within 1270 — 1480 A. This division allowed to see weak emission lines, which 
were not detected in previous studies of the quasar composite spectra, but were 
known from individual optical or composite UV spectra from the Hubble Space 
Telescope. Although the physical explanation of the difference in spectral indices 
between quasars and their dependence on quasar parameters is still not clear, it is 
obvious that this difference has to be taken into account when generating composite 
spectra, e. g. for redshift measurements. It was also shown that the equivalent 
width of the emission lines does not depend on the spectral index. 



Introduction 



The spectral energy distribution of quasars in UV-optical range is characterised by the so called 
Big Blue Bump with a pick around X res t — 1000 — 1300 A, broad emission lines of hydrogen Lyman 
and Balmer series and low-ionized metals (e.g. Civ, Sim, Mgn etc.), broad absorption lines (only in 
~ 10% of quasars) and decrement of the flux blueward of 1215 A due to absorption by intergalactic Hi 
(the Lycv forest). It is believed that the (mostly thermal) continuum and emission lines originate from 
the hot accretion disk and circumnuclear fast moving clumps of gas, correspondingly. The proximity o f 
these regions is considered to be the most promising explanation of the Baldwin effect Baldwin! . Il977j | : 
the inverse correlation of equivalent width of some emission lines with the monochromat ic fluxes at UV 



region (in most cases the flux at 1450 A or 2500 A is considered in the literature [e.g. iDietrich et al 



20021 : IWu et all . l2009j 



Due to remarkable similarity of the quasar spectra from object to object, the mean, or composite, 
spectra are usually used to study the general spectral properties of quasars. They are also widely 
used as templates for redshift determination (mainly in automatic redshift surveys). The averaging 
over several hundreds or even thousands of (usually) medium-resolution spectra allows to increase the 
signal-to-noise ratio and to reveal tiny features unseen in individual spectra. The co mposite spectra of 
quasars in ultraviolet-opt i cal band were compiled for a wide set of qu asar samples [Brotherton et al 



20011 : iFrancis et all Il99lt IScott et all 12004 Ivanden Berk et all l200ll. Thev also used i n Lya-forest 
studies [e.g. Bernardi et al.l ~ 2003 : Desiacaues et all 2007 ; Polinovskvi and Malygin . 201d| . One of the 
main sources of uncertainties during the generation of the quasar composite spectra is the difference in 
spectral shape of individual objects (despite of the general similarity), e.g. the difference in spectral 
index a\ of the red part of UV- optical bump. Thus in the present work to study the part of the 
UV-optical bump we used the composite spectra compiled by us from samples of spectra with similar 
ax- 



The data 



For our study we used 16 composite spectra of quasars compiled by us Ivashchenko et all 120121 ] . 



Each of these spectra was generated of 200 individual ones with similar spectral indices a\ within the 
range 1270 — 1480 A namely from those that have the values of ct\ close to —0.6 — k ■ 0.1, where k — 1..16 
is the composite number. T hese subsamples were ta ken from the publicly available release of the sky- 
residual subtracted spectra Wild and Hewett . 2010l j for the Sloan Digital Sky Survey (SDSS) Legacy 
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Release 7 [Abazaiian et al. , 2009j . The spectral indices for individual spectra were calculated considering 
the continuum of the quasar spectra redward of the Lya emission line (1215 A) to be a power-law ~ A" A , 
and using the following wavelength ranges, which are the most free from emission lines: 1278—1286, 
1320-1326, 1345-1360 and 1440-1480 AA. The mean arithmetic spectra were obtained after (i) division 
of each individual spectrum on its normalization constant (the mean flux over all pixels within the rest 
wavelength range 1450-1470 A), (ii) reb inning them with AA res t = 2 A, and (hi) stacking the spectra 
into the rest frame. The dispersion a 1 of each pixel of the composite spectrum is calculated from the 
noises o~i of pixels of individual spectra as er~ 2 = ^2 • 

i 

The spectral indices of the composite spectra were calculated in the way mentioned above, but the 
parts of the spectra which are the most free from emission lines were selected manually in each composite. 
The parts of spectra used for calculating the spectral index are shown in Figure [1] (short-dashed line, 
the highest spectra are the steepest ones). Three spectra with the slopes —0.84, —1.49 and —2.12 are 
shown with the solid lines together with the fitted power-law continua (dashed lines). The regions used 
for continuum fitting slightly vary from spectrum to spectrum which is probably caused by changing of 
intensity of weak emission lines. It is worth to note that there is no absolutely emission line-free region in 
UV-optical bump, and selection of lines-free regions means selection of regions with the weakest emission 
lines. 
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Figure 1. The parts of the composite spectra used for continuum approximation. See explanation in 
the text. 



Table 1. The mean redshift z, spectral index a\, and the mean monochromatic luminosity (log/1450) 
at 1450 A for 16 composites. 
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The spectral indices for the composite spectra, the mean redshift of the corresponding sample, and 
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the mean logarithm of the quasar monochromatic luminosity at 1450 A, (logZi45o), are presented in 
Table 1. The the value logZi45o for individual spectra was calculated using the mean flux within the 
wavelength range 1449 — 1451 A and photometric distance to the quasar obtained within the spatially 
flat ACDM cosmological model with the matter density parameter fijv/ = 0.27 and Hubble constant 
Hq = 70.5 km s -1 Mpc -1 . The quoted error bars are the root mean squares for the corresponding 
distributions for z and (log^so), and la errors for a\, correspondingly. 



The method 

The wavelength ranges (a) 1215-1321A and (b) 1323-1447A (Figure© were considered separately 
and fitted with a superposition of continuum and the smallest possible number of emission lines in the 
form: 

(A'' est - A°) 2 



/(A rest ) = C + ^a iex p 



2wt 



(1) 



where in case (a) C = b, i. e. the continuum is fitted with a constant, and in case (b) C = d\ ax with 
fixed a\ from Table 1. The values of b, d, a,i, Wi, A° were calculated in the following two steps: (i) 
using the IDL lmf it subroutine for each multiplet the best fit model (based on the minimal x 2 /d.o.f.) 
in a form of ((J) was found; (ii) the central wavelengths A obtained at the first step were fixed and 
the best fit values of other parameters {b,d,ak,Wk} with the la errors from the extremal values of the 
N-dimensional di stribution were ca l culate d by the Markov Chain Monte Carlo (MCMC) method using 
CosmoMC package Lewis and Bridld . l2002j | as a generic sampler (with the values of {b, d, a^., Wk} obtained 



at the first step as starting values). The MCMC technique has been chosen as it is a fast and accurate 
method of exploration of high-dimensional parameter spaces. In each case we generated 8 chains which 
have converged to R — 1 < 0.0075. Due to the small values of the flux dispersion a"j in composite spectra 
determined from the covariance matrix to estimate errors correctly we introduced the additional intrinsic 
errors a int such that the total dispersion, a 2 = a 2 - + af nt , results in the minimal x 2 /d.o.f. to be 1 — 0. 

For separate lines or sum of lines the equivalent widths were calculated. In each case the area under 
the lines was calculated as an integral of an analytic function (Gaussian or several Gaussians) with the 
obtained parameters. 



Results and discussion 

In Figure [3] one of composite spectra with a\ = —0.92 is shown with the best fits for both regions 
as an example. The dashed lines indicate the emission lines found in spectrum. The labels X1-X5 
stand for components, identification of which was not found in the literature. In Table [2] the central 
wavelengths of emission lines from all 16 spectra are presented along with possible identification and cor- 
responding laboratory rest wavelengths. The obtained parameters for each line (amplitude and FWHM) 
will be presented elsewhere. The identificatio n of lines was made using known emission lines from the 
literature. The following lines were found by IVanden Berk et al. 2001] in composite spectrum of 2200 



quasars from SDSS without taking into account any differences in spectra shape: Lya+Oiv (one fea- 
ture), Nv, Sin, Oi+Sin (one feature), Cn and Siiv+Oiv] (one feature). The separate components of 
these emiss ion features as far as other lin es found in the present work are known f rom works on ultraviolet 
individual |Laor et aL , Il994l 1995 , Il997 1 and composite [Brotherton et all Il994 1 spectra of quasars from 



the Hubble Space Teles c ope, and from empirical UV template for Fe emission in quasars obtained by 
Vestergaard and Wilked 2001 1 from the galaxy I Zwicky 1. All the lines, for which we did not find any 



possible candidates in the literature could be new lines, not noticed by other authors due to their small 
intensity, as well as additional components of the known lines originated due to some special physical 
conditions in the broad line regions in a quasar. 

In Figure [4] the spectral index - equivalent width (for separate lines or sum of lines) diagrams are 
presented. We applied the F-test and found that for all emission features presented in this figure the 
values of equivalent width do not depend on the spectral index. This result is an expected one if we 
take into account the Baldwin effect and the absence of luminosity-dependence of the spectral index, as 
it is seen from Table [TJ On the other hand, the mean redshift of the subsamples slightly increases with 
a\, that could be an evidence for evolution of a\ with time, but also can be caused by some selection 
effects. Although the physical explanation of the difference in spectral indices between quasars and their 
dependence on quasar parameters is still not clear, it is obvious that this difference has to be taken into 
account when generating the composite spectra, e. g. for redshift determination. When using composite 
spectra made without taking into account the difference in spectral shape one usually deals with larger 
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Figure 2. Wavelength ranges 1215-1320 A (left) and 1323-1443 A (right). Dashed lines indicate the 
central wavelengths of the emission lines found by iVanden Berk et ail 2001j . Here asterisk indicates 
more than one allowed value of the total angular momentum J for that specific term and transition, "]" 
means intcrcombination transitions. 



noise and thus treats blended emission features as one line, increasing the errors in central wavelengths 
and, hence, in line identification. 




1250 1300 1350 1400 1450 

A, angstroems 



Figure 3. The composite spectrum with a\ = —0.92 with the best fits for both two regions. The dashed 
lines indicate the found emission lines, while short-dashed lines indicate their central wavelengths. XI - 
X5 stand for lines which were not identified by us. 
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Table 2. The cental wavelength of emission lines found in spectra. Lines in brackets are probably 
components of one found emission feature. XI etc. stand for lines which were not identified by us. Here 
asterisk indicates more than one allowed value of the total angular momentum J for that specific term 
and transition, "1" means intercombination transitions. 
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